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SUMMARY 

DTIC  QUALITY  mSPECTED  3 

Bonded  fibre -composite  doublers  are  increasingly  being  used  to  reinforce  andlor  repair 
damaged  or  underdesigned  metallic  aircraft  structure.  This  trend  will  continue  as  the 
average  age  of  Australian  civilian  and  military  aircraft  increases.  Two  mathematical  models 
of  the  stress  state  in  bonded  doublers  are  presented  in  tins  report:  the  end-tapered  double-lap 
joint  and  the  stepped  double-lap  joint.  The  mathematical  development  of  the  two  models  is 
detailed  and  the  fidelity  of  the  predicted  stress  states  compared  with  that  obtained  from  finite 
element  analyses  and  experiment. 
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1.  INTRODUCTION 

Aging  aircraft,  both  civilian  and  military,  are  prone  to  various  tj^es  of  deterioration  and 
damage  in  their  metallic  components  including  general  corrosion  and  cracking  due  to  stress 
corrosion  and  fatigue.  A  technique  pioneered  at  ARL‘  to  extend  the  life  of  these  aircraft  is 
to  bond  fibre-composite  reinforcements,  such  as  boron/epoxy  or  graphite  epoxy,  over  the 
damaged  regions.  In  demanding  situations,  such  as  highly-stressed  thick-section  repairs, 
accurate  methods  of  design  and  analysis  are  needed.  This  report  details  the  development  of 
two  mathematical  models  for  the  adhesive  shear-strain  distribution  in  a  scarfed  repair;  the 
end-tapered  double-lap  joint  (Figure  la)  and  the  stepped  double-lap  joint  (Figure  lb). 
These  two  joints  are  idealised  representations  of  an  end-tapered  repair.  An  actual 
end-tapered  repair  or  doubler,  such  as  the  ARL^  developed  F-1 1 1  wing-pivot-fitting  (WPF) 
doubler  (Figure  2),  has  a  geometiy  as  shown  in  Figure  Ic,  where  after  lay-up  the  doubler 
has  been  inverted  and  the  outer  plies  bent  over  by  pressure  in  an  autoclave.  Also  in  this 
report  is  a  comparison  between  the  mathematical  model  developed  herein  and  a  finite 
element  (FE)  analysis  for  both  the  end-tapered  and  the  stepped  double-lap  joint.  Finally  a 
comparison  is  made  between  the  measured  strains  in  the  outer  adherend  of  a  boron  doubler 
(Figure  Ic)  and  the  strains  predicted  by  the  models  of  the  end-tapered  and  the  stepped 
double-lap  joints. 


(a)  cnd-tapered  double-lap  joint 


=1 


(b)  stepped  double-lap  joint 


(c)  schematic  of  actual  joint 
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The  primary  reason  for  tapering  the 
ends  of  a  composite  doubler  is  to 
reduce  the  peel  and  shear  stress,  which 
arises  there  due  to  the  load  transfer 
conditions.  A  high  peel  stress  (tensile 
stress  that  acts  transverse  to  the  plane 
of  the  bond)  is  detrimental  because  the 
interlaminar  tensile  strength  of  a  fibre 
composite  lay-up  is  typically  less  than 
10%  of  its  longitudinal  tensile  strength. 

For  example,  the  longitudinal  tensile 
strength’  of  unidirectional  boron/epoxy 
is  1260  MPa  whereas  its  interlaminar 
tensile  strength  is  only  61  MPa  -  and 
the  ratio  of  interlaminar  to  longitudinal 
toughness  is  even  smaller  than  that  of 
the  strengths.  Adhesives  are  also 
susceptible  to  failure  induced  by  peel 
stresses.  Tapering  the  ends  reduces  the 
adhesive  shear-stress  concentration  at 
the  ends  of  the  doubler  and  hence  also 
reduces  the  peel  stress  resulting  from 
the  unbalanced  shear  stress  acting  on 
the  doubler.  An  optimum  taper  angle  is 

generally  about  three  degrees  since  this^®^®******®  Two  Mathematical  Models  and 
gives  an  adequate  degree  of  stress  relief  Actual  Doubler 

while  ensuring  the  size  of  the  doubler  is 
not  too  great.  The  modelling  of  the 
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aluminium  inner  adherend 
adhesive 

fibre-composite  outer  adherend 


Figure  1 


2 


peel-stress  distribution  is  a  more  difficult  mathematical  problem  and  is  not  addressed  in  this 
report,  however,  minimising  the  shear  will  also  minimise  the  peel  stress. 


The  application  of  doublers  to  an  aircraft  structure  can  induce  some  bending  due  to  the 
added  eccentricity  of  the  load  path.  However,  most  aircraft  structures  have  supports  that 
act  to  reduce  this  bending.  For  example,  the  WPF  of  the  F-1 11  (Figure  2)  has  stiffeners 
and  shear  webs  and  is  bolted  to  the  sub-structure.  The  strength  of  the  single-sided 


OUTBOARD  < 


Figure  2.  F-1 11  Doubler-Section  Looking  Aft  on  Right  Wing. 


end-tapered  joint  will  then  approach  that  of  the  end-tapered  double-lap  joint.  Analysis  of 
the  single-sided  end-tapered  double'  can  then  be  approximated  by  considering  either  the 
idealised  geometry  of  the  end-tapered  double-lap  joint  (Figure  la)  or  the  stepped 
double-lap  joint  (Figure  lb)  and  by  using  a  'mechanics  of  solios'  approach  similar  to  that 
employed  by  Hart-Smith''  in  the  analysis  of  the  double-lap  joint 

2.  THE  END-TAPERED  DOUBLE-LAP  JOINT 

2.1.  Elastic  Analysis 

2.1.1.  Analytical  Development 


In  this  analysis,  it  is  assumed  that  the 
adhesive  deforms  only  by  shear  and  that  the 
adherends  deform  only  by  stretching  in  the 
longitudinal  direction.  This  one-dimensional 
approach  works  reasonably  well  for  the 
non-tapered  double-lap  joint^  and  should, 
given  the  reduced  peel  stresses,  work  even 
better  for  the  end-tapered  double-lap  joint. 
Geometry  and  nomenclature  for  the 
end-tapered  double-lap  joint  is  shown  in 
Figure  3.  A  fibre-composite  outer  adherend 
can  be  accommodated  by  using  an  effective 
modulus  in  the  longitudinal  direction.  As 
most  doublers  employed  at  ARL  have 
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Figure  3.  Geometry  and  Nomenclature 
for  the  End-tapered  Double-Lap  Joint 
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unidirectional  reinforcement,  the  modulus  along  the  fibre  direction  of  a  single  ply  can  be 
used.  To  aid  the  analysis,  the  joint  is  divided  into  two  regions:  the  tapered  region  from  jc=0 
to  J£:=l  and  the  non-tapered  region  from  x=\  to  (or  4  =0  to  4  =«)•  The  differential 
equation  describing  adhesive  deformation  in  the  first  region  is  derived  below  and  the 
differential  equation  for  the  second  region  can  be  found  from  the  first  by  putting  the  angle 
of  taper  to  zero. 


Force  equilibrium  in  the  horizontal  direction  dictates  that; 


£[0 

dx 


-x  =  0. 


(1) 


and, 


£Ti 

dx 


+  2x  =  0. 


(2) 


Here,  x  is  the  adhesive  shear  stress  and  T,  and  T,  are  tractions  in  the  outer  and 
inner  adherends  respectively.  The  elastic  stress-strain  relations  for  the 
adherends  are. 


dx 


To 


Eo^jc  tan0  -4-  to° 


(3) 


and, 


dx  Eitj' 


(4) 


The  terms  5„  and  5,  are  the  displacements  of  the  outer  and  inner  adherends 
respectively.  The  elastic  shear-stress  shear-strain  relation  for  the  adhesive  is. 

X  =  Gy  .  (5) 

In  this  analysis,  the  adhesive  shear  strain  is  averaged  across  the  thickness  of  the 
adhesive  layer; 


r= 


So  -  Si 

r) 


(6) 


Taking  the  derivative  of  equation  6  with  respect  to  (w.r.t.)  jc  and  substituting  in 
equations  3  and  4  yields. 
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dy 

dx 


1 

T1 


To 


Eo  X  tan6  + 


(7) 


Differentiating  again  w.r  t.  x  yields. 


d^y  ^  1 

dx'^  ^ 


Eo^x  tan0+  to^^- EoTotanB 
Eo^^jc  tan0+  to^ 


£EiJ_ 

dx  Ejti 


(8) 


but 


equations  1  and  2  imply  T„  =  X  dx  and  — — ^  =  -2x,  so  that  equation 

•'0  dU: 


8  becomes  in  terms  of  shear  stress 


d 

dx'^ 


1 

n 


Eo^j:  tan0+  to^T-  EotanO j” t(jc )tix 

to») 


l' 

Eo^f  tan0+  tc 


+  2X: 


.  (9) 


Multiplying  both  sides  of  equation  9  by  Eo^|^JC  tan0+  J  ,  substituting  X  =  G^, 
and  taking  the  derivative  w.r.t.  x  yields  the  third-order  ordinary-differential  equation: 


d^y  ^  2tan0  d 

dx^  jc  tan0+  to  dx^ 


G 

T1 


Ei  ti 


+ 


1 


Eo  X  tan0+  to 


d  y 
dx 


4tan0 


Ei  ti^x  tanO-b  to  j 


Y- 


(10) 


Equation  10  is  probably  not  soluble  in  closed  form.  It  is  amenable,  however,  to  a  numerical 
solution  using  the  NAG'  subroutine  D02HAF.  This  is  a  subroutine  for  solving  boundary-value 
problems  involving  ordinary  differential  equations  of  arbitrary  order.  The  equations  are 
re-written  as  a  system  of  simultaneous  first-order  differential  equations  and  the  problem  solved 
by  a  shooting  and  matching  technique. 
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Numerical  Algorithms  Group  library  of  mathematical  subroutines. 
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The  boundary  condition  at  jc=0  is  obtained  from  equation  7; 


P 

TlEi  ti‘ 


(11) 


The  other  two  boundary  conditions  are  derived  by  matching  y  and  y*  at  x=l  from 
the  preceding  solution  with  y  and  y*  at  1^=0  from  the  solution  for  the 
non-tapered  double-lap  joint.  The  form  of  the  solution  for  the  non-tapered 
double-lap  joint  is: 


(12) 


where  the  subscript  n  refers  to  the  non-tapered  double-lap  joint  and  P,  is  the 
load  remaining  at  jc=I  to  be  transferred  from  the  inner  to  the  outer  adherends  at 
JC=I. 

The  solution  shown  above  for  the  non-tapered  double-lap  joint  can  be  obtained 
from  the  previous  equations  by  substituting  6  =  0  and  in  equation  8. 

Then. 


_  1 


1  dTo  1  dTi 


dx'^  'HLEoto  dx  Ejt,  dx  j' 
Substituting  the  equilibrium  equations  1  and  Z  and  equation  5  gives; 


(13) 


d^yn  ^  G 
dx^  n 


LEoto°  E.t,  J 


Yn, 


=  P^Yn. 


(14) 


whefp 


(15) 


The  solution  of  equation  14  over  a  semi-infinite  domain  is  equation  12  above. 

The  boundary  condition  y(«>o)  =  0  is  implicit  in  this  formulation.  Solution  over  a 
finite  domain,  p  say,  is. 

Yn  =Acosh  (^p(^p  -  ^  j  ^Bsinh  (^p^p  -  ^  j  (16) 

P,  in  equation  12  is  the  load  remaining  at  jp'I  to  be  transferred  from  the  inner  to 
the  outer  adherends  and  is  given  by: 
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Pi=Ti-P 


Ej  tj 

2Eo  to  -t-  Ei  ti 


(17) 


The  last  term  in  equation  17  expresses  the  condition  that  at  jc  =  «  the  load  is 
distributed  to  each  adherend  according  to  its  stiffness. 

Matching  at  j:=I  (4=0)  implies  that: 


Y(I)  =  Yn(0). 


and. 


Pi 

TlE.t.p’ 


(18) 


P| 

nE.t, 


(19) 


Since  P,  is  not  known  a  priori  neither  are  the  boundary  conditions  at  jc=I. 
Consequently,  an  iterative  scheme  was  implemented  on  computer  to  solve  the 
governing  differential  equation  (equation  10)  and  its  boundary  conditions 
(equations  11,  18,  and  19).  The  iterative  scheme  involved  making  repeated 


guesses  at  the  parameter  P,  until  the  load  transferred  to  the  outer  adherends 

Eo  t 

(which  was  found  by  integrating  equation  1)  was  equal  to  P  — . 

to  "h  fc;  tj 


2.1.2.  Comparison  of  the  Elastic  Model  with  an  FE  Analysis 


An  FE  analysis  of  an  end-tapered  double-lap  joint  was  set  up  in  PAFEC.  Only  one  quarter 
of  the  joint  needed  to  be  modelled  due  to  the  twofold  symmetry  of  the  joint.  The  adhesive 
and  the  adherends  were  modelled  using  eight-noded  isoparametric  elements.  Previous 
experience  in  modelling  double-lap  joints  dictated  the  need  for  four  elements  through  the 
thickness  of  the  adhesive  layer  to  achieve  convergence.  Adhesive  shear  strains  were 
calculated  by  subtracting  the  longitudinal  displacements  of  nodes  on  either  side  of  the 
adhesive  layer  and  dividing  through  by  the  adhesive  thickness.  The  adhesive  shear  strain  is 
therefore  averaged  across  the  thickness,  as  in  the  analytical  development,  thus  allowing 
proper  comparison.  The  graphs  that  follow  show  the  effect  of  varying  joint  parameters  on 
the  adhesive  shear-strain  distribution  in  the  joint.  Table  1  contains  those  joint  parameters. 
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Figure  7. 

Figures  4,5,6,  and  7  show  that  the  model  of  the  end-tapered  double-lap  joint,  assuming 
elastic  adhesive  deformation,  is  in  excellent  agreement  with  the  finite-element  analysis  for  a 
wide  range  of  joint  parameters. 

2.1.3.  An  Approximate  Analytical  Solution:  a  First-Order  Perturbation  Solution 


Distance  along  Joint  (mm) 
Figure  6. 


.\n  approximate  analytical  solution  to  the  differential  equation  governing  elastic 
deformation  of  the  adhesive  (equation  10)  can  be  found  by  proposing  a  regular  perturbation 
expansion  of  the  shear  strain; 
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7(^)  =  yoix)  +  eyi  (x) + e^y2(x)  + . 


where  the  parameter  e  is  a  small  number  much  less  than  one.  This  parameter  is 
associated  with  term  tanB  in  equation  10.  If  the  taper  angle,  0  ,  is  very  small 
then  tanQ  will  also  be  small.  Rewriting  equation  10  in  terms  of  e  gives; 

The  first  order  perturbation  expansion  of  the  shear  strain  (i.e.  the  first  two  terms 
on  the  right  hand  side  of  equation  20)  is  then  substituted  into  equation  21  and 
terms  with  like  powers  of  e  collected.  The  ztroth  order  terms  are; 

The  relevant  boundary  conditions  are; 

^(0)  =  ^ 
dx  tiEit,’ 


The  parameter  P  ,is  given  by: 

The  zeroth  order  function,  y„,  is  then: 

Yo(x)  =A  cosh^Po^I-xjj+Bsinh^Po^I-xjj^ 

The  parameter  P  „  is  given  by  . 


The  constants  A,  B,  C  are  found  by  substituting  the  boundary  conditions  into 
equation  23 


llEitiPo’ 


Note  that  the  parameter  P,  is  not  yet  determined.  It  is  found  by  setting  the 
integral  of  the  shear  stress  equal  to  the  load  per-unit-width  that  has  to  be 
transferred  from  the  inner  to  each  outer  adherend. 

The  first  order  terms  found  on  substituting  the  perturbation  expansion  into 
equation  10  are: 


The  relevant  boundary  conditions  are; 


7,(l)  =  0; 

y{(0)  =  0; 

The  homogeneous  solution  to  the  above  differential  equation  is; 

y;'^(x)  =Dcosh[^Po^l  -X  j  j+Esinh[^po(^l  -x)]+F.  (25) 

The  particular  solution  to  the  above  differential  equation,  having  first 
substituted  for  Yo(x),  is; 

=aa  x+bb  X  cosh  [p„(l  -  x]]kc  x  sinh  [po(l  -  x]]4dd  x*  cosh[po(l  -  x^] 
+CC  x^  sinh^Po^l  -X  j  j.  (26) 


The  constants  aa,  bb,  e.t.c.  are  given  by; 
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Figure  8  shows  a  comparison  of  the  numerical  solution  for  the  joint  configuration  that 
follows  and  the  first-order  perturbation  solution. 


12 


distance  along  bondline  normalised  against  taper  length 


Figure  8. 


T1 

(mm) 

G 

(MPa) 

Eo 

(MPa) 

c 

(mm) 

E. 

(MPa) 

t, 

(nun) 

e 

(degrees) 

1 

(mm) 

P 

(N/mm) 

0.14 

700 

208e3 

0.13 

71e3 

6.36 

2.42 

27 

500 

Table  2. 

The  1st  order  perturbation  solution  shows  reasonable  agreement  with  the  more  accurate 
numerical  solution.  A  higher  order  perturbation  expansion  for  the  adhesive  shear  strain 
would  be  necessary  to  achieve  a  more  accurate  approximation. 

2.2.  Elastic/Perfectly-Plastic  Analysis 
2.2.1.  Analytical  Development 

To  obtain  the  solution  for  the  adhesive  shear-strain  distribution  in  a  joint  deforming  in  an 
elastic/perfectly-plastic  manner  the  tapered  region  must  be  divided  into  an  elastic  region 
and  a  plastic  region.  Firstly  the  adhesive  shear-stress  shear-strain  curve  is  idealised  as  being 
elastic/perfectly-plastic. 

In  the  plastic  region  a  solution  is  obtained  by  setting  T(x)»Tp(the  yield  stress)  in 
equation  9: 
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£l=i. 

dx'^  ^ 


T  E  t° 

_ Wp  1^0  ip 

tan8  +  to 


Eiti 


Integrating  equation  30  w.r.t.  x  yields; 


^_1 

dx  ^  Ei  ti 


Eotan6  jctan0  +  to 


Integrating  once  again: 


y{x)  =  ^  ^ 
T1  Eiti 


Xp  t“ln^jctan0  +  t^^ 

Eotan20 


+  Ok  +  D. 


Applying  the  boundary  condition  given  by  equation  11  to  equation  32  yields; 


T1  Eiti 


2  Xpto  ln^xtan0  +  to  j 


Eo  tan^0 


- E-L  +  d 

T]Eotan0  TlEitiJ 


As  in  the  elastic  analysis  the  other  boundary  conditions  are  found  by  matching.  In  this  case, 
however,  the  adhesive  in  the  tapered  region  can  deform  fully  plastically  or  partly  plastically 
and  partly  elastically.  If  the  deformation  in  the  taper  is  fully  plastic  then  equation  33  is 
matched  with  the  solution  describing  adhesive  plasticity  in  the  non-tapered  double-lap  joint 
at  x=f.  That  solution  is  obtained  by  setting  tan6  =  0  in  equation  30  and  integrating; 


fe-lLr_L_  +  2 

d^2“nLEoto  EitJ  G' 


and  hence. 


-lEg!V 


hE^+F. 
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This  solution  is  then  matched  with  that  assuming  adhesive  elasticity  in  the  non-tapered 
double-lap  joint  (equation  12)  at  the  plastic-to-elastic  transition  point.  The  location  of  this 
transition  point  is  not  known  a  priori  and  so  an  iterative  scheme  was  again  implemented  on 
computer  to  find  the  full  solution  by  making  guesses  at  the  starting  strain. 

If  the  adhesive  deforms  plastically  for  only  part  of  the  scarfed  region  then  a  different 
procedure  is  necessary.  Equation  33  is  now  matched  with  the  equation  10  (the  solution 
assuming  adhesive  elasticity  in  the  scarfed  re^on)  at  the  plastic-to-elastic  transition  point. 
At  x=l  equation  10  is  matched  with  the  solution  assuming  adhesive  elasticity  in  the 
non-tapered  joint  (equation  12).  A  computerised  iterative  solution  was  again  necessary  to 
find  the  full  solution. 

2.3.  An  Upper  and  a  Lower  Bound  for  the  Peak  Shear  Strain 

It  is  desirable,  from  an  engineering  design  viewpoint,  to  find  a  simple  expression  for  both 
the  upper  and  lower  bound  for  the  peak  shear-strain  in  an  end-tapered  double-lap  joint.  A 
simple  approach  is  to  consider  two  separate  double-lap  joints.  The  joint  used  to  estimate 
the  upper  bound  would  have  an  outer  adherend  thickness  equal  to  t.  and  the  joint  used  to 
estimate  the  lower  bound  would  have  an  outer  thickness  of  t,®.  By  comparing  the  results 
from  the  end-tapered  double-lap  joint  model  (with  fixed  values  of  and  t,®  but  varying 
values  of  0)  with  the  results  fi'om  the  two  double-lap  joints  the  usefulness  of  the  estimates 
for  the  upper  and  lower  bounds  can  be  evaluated.  To  be  useful,  at  small  values  of  6  the 
peak  shear-strain  should  approach  the  lower  bound  and  at  large  values  of  0  it  should 
approach  the  upper  bound.  Two  cases  were  investigated;  elastic  and 
elastic/perfectly-plastic  adhesive  deformation. 

2.3. 1 .  Assuming  Elastic  Adhesive  Deformation. 

Figure  9  shows  the  variation  of  peak  shear  strain  with  taper  angle  when  the  adhesive  is 
assumed  to  deform  elastically.  The  joint  parameters  are  as  for  Figure  4  with  these 
exceptions:  P=2000  N/mm  and  Ti=0. 14mm. 


Figure  9. 
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The  upper  bound  provides  a  sound  upper  limit  for  the  peak  shear  strain  at  large  angles  of 
taper  but  the  lower  bound  underestimates  the  peak  shear  strain  at  small  angles  of  taper. 
Taper  angles  in  a  repair  are  typically  2-4°  and  hence  the  upper  bound  may  be  too 
conservative  to  be  used  for  design  purposes.  Likewise,  the  lower  bound  may  be  too  low  to 
be  used  for  design  purposes  (e.g.  for  a  taper  angle  of  2.7°  the  actual  peak  shear  strain  is 
58%  greater  than  the  lower  bound). 

2.3.2.  Assuming  Elastic/Perfectly-PIastic  Adhesive  Behaviour 

When  the  adhesive  is  allowed  to  deform  in  an  elastic/plastic  manner  a  similar  graph  is 
obtained  (Figure  10).  In  this  case  the  lower  bound  greatly  underestimates  the  peak 
shear  strain.  Loads  high  enough  to  cause  adhesive  plasticity,  however,  are  rarely 
encountered.  Table  3  contains  the  parameters  for  the  joint. 


p 

(N/mm) 

Ei 

(MPa) 

ti 

(mm) 

Eo 

(MPa) 

C 

(mm) 

to 

(mm) 

G 

(MPa) 

n 

(mm) 

4,000 

71,000 

6.36 

205,000 

0.13 

1.27 

750 

37 

0.13 

Tables. 
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Figure  10. 

3.  THE  STEPPED  DOUBLE-LAP  JOINT 

3.1.  Elastic  Analysis 

3.1.1.  Analytical  Development 

This  model  of  the  stepped  double-lap  joint  (Figure  lb)  applies  the  solution  for  the 
double-lap  joint  (equation  1 6)  to  each  step.  Each  ply  of  the  outer  adherend  is  assumed  to 
be  of  step  length  I,  thickness  t,,  and  Young's  modulus  E,.  A  joint  of  this  type  could  be 
made  by  stepping  the  plies  of  a  fibre  composite  lay-up  to  form  the  outer  adherends.  The 
form  of  the  solution  on  the  ^h  step  is. 
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Yk  (4)  =AkSinh  j+BtCosh 


where  ^  is  the  local  coordinate  for  each  step  and  where 

The  boundary  condition  given  by  equation  11  applies  to  the  first  step; 

^(0)  =  — P_ 

BEiti 


=A,P, 


At  the  end  of  each  step,  of  length  I,  continuity  of  the  shear  strain  must  be 
enforced,  though  the  derivative  of  the  shear  strain  may  not  necessarily  be 
continuous  since  the  profile  of  the  outer  adherend  is  not  smooth. 


Therefore, 


Yk(l)  =  Yk+i(0) 


Bk+i=AkSinh  J^Pkl  j+BkCOsh  (^Pklj. 


At  the  start  of  each  step,  the  derivative  of  the  shear  strain  is: 

^(0)  =  ir— 12 - 

dC  ^Uotok  Eitij- 


This  equation  can  be  derived  from  equation  7  by  setting  tanS  =  0  and  tS=kto  • 

If  the  load  transferred  from  the  inner  to  each  outer  adherend  over  the  kth  step  is 
then  at  the  start  of  the  k+Ith  step  the  tractions  are: 


T„=2;  Fj, 

i=i 
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Ti=P-2£Fj . 

j=i 


Therefore,  from  equation  41; 


AkPk 


dy,  i(  IF3  P-ISF, 


"nvEotok  Eiti 


where  F,^  is  given  by  integrating  equation  36: 


Fk  =G 


Akj^cosh  ^pklj-  Ij  +  BkSinh^Pkl 


(43) 


(44) 


(45) 


Thus  if  B,  in  equation  36  is  known  then  ail  the  coefficients  and  can  be 
found.  At  the  start  of  the  last  step  the  shear  strain  must  also  be  continuous. 
The  form  of  the  solution  for  this  last  («th)  step  is  as  for  equation  12: 


(46) 


Continuity  of  the  shear  strain  implies; 


tie 


=An.l  sinh  1  y  B „.1  cosh  (^pn-l  1  j • 


(47) 


The  total  load  transferred  from  the  inner  adherend  to  each  of  the  outer 
adherends  is: 


2:f,+ 


/=! 


Pi 

EitiPn'' 


(48) 


Therefore,  all  that  remains  is  to  find  B,  and  then  all  other  coefficients  can  be  found  from 
equations  41  and  45  and  P,  can  be  found  from  equation  47.  A  computerised  iterative 
scheme  was  again  implemented  in  which  repeated  estimates  were  made  of  Bj  and  the  load 
transferred  to  the  outer  adherends  calculated  (from  equation  48)  until  that  value  was  to 


equal  to  P 


_ Eq  tp 

2Eon  to+Ejtj  ■ 


3.1.2.  Comparison  of  the  Elastic  Model  with  an  F£  Analysis 


An  FE  analysis  of  a  stepped  double-lap  joint  was  set  up  in  PAFEC.  The  model  was  similar, 
in  terms  of  element  type  and  spacing,  to  the  end-tapered  double-lap  Joint  model.  Once 
again,  only  one  quarter  of  the  joint  n^ded  to  be  modelled  due  to  the  symmetry  of  the  joint. 
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The  adhesive  shear  strain  was  again  calculated  by  subtracting  displacements  across  the 
adhesive  layer  and  dividing  by  the  adhesive  thickness.  Three  joints  were  modelled  and  the 
table  of  their  joint  parameters  and  figure  numbers  are  shown  below.  The  figures  show 
excellent  agreement  between  the  mathematical  model  and  the  FE  results. 


Figure  1 1 

Figure  12 

Figure  13 

E„(MPa) 

208,000 

208,000 

138,000 

t„  (mm) 

0.13 

0.13 

0.15 

E.  (MPa) 

71,000 

71,000 

71,000 

t.  (mm) 

6.36 

6.36 

6.36 

n  (no.  steps) 

10 

10 

10 

1  (mm) 

3 

5 

5 

P  (N/mm) 

2,000 

2,000 

2,667 

G(MPa) 

590 

590 

590 

T)  (mm) 

0.1 

0.2 

0.2 

Table  4. 


Figure  1 1 . 


distance  along  joint  (mm) 


Figure  12. 
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distance  along  joint  (mm) 
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3.1.3.  An  Upper  and  a  Lower  Bound  for  the  Peak  Shear  Strain 

Bounds  for  the  peak  shear  strain  in  the  adhesive  assuming  elastic  deformation  ca.a  be 
developed  in  much  the  same  way  as  for  the  end-tapered  double-lap  joint.  The  lower  bound 
is  taken  to  be  the  adhesive  shear  strain  when  only  one  outer-adherend  ply  is  present.  The 
upper  bound  is  taken  to  be  the  adhesive  shear  strain  when  all  plies  are  present  in  a 
non-stepped  geometry.  Table  five  below  has  the  parameters  for  this  joint.  The  results  are 
shown  in  figure  14. 


P  (N/mm) 

Ei  (MPa) 

ti  (mm) 

Eo  (MPa) 

to  (mm) 

n 

G(MPa) 

Ti  (mm) 

1,500 

71,000 

6.36 

205,000 

0.13 

5 

750 

0.13 

Table  5. 


The  lower  bound  is  quite  sound  for  this  joint  configuration.  The  peak  shear  strain  for  this 
joint  with  a  step  length  greater  than  6  mm  can  be  readily  approximated  by  a  joint  having  a 
single  outer-adherend  ply.  This  step  length  is  equivalent  to  a  taper  angle  of  about  1  2°.  A 
typical  step  length  is  about  3  mm  (2.6°)  in  which  case  the  lower  bound  underestimates  the 
peak  shear  strain  by  about  20%. 


0  12 
011 
010 

peak 
shear 

strain  oo? 

000 
006 
OM 

0  2  4  €  <  fO  f4 

Step  length  (mm) 

Figure  14. 

3.2.  Elastic/Perfectly-Plastic  Analysis 
3.2.1.  Analytical  Development 

Solving  for  the  elastic  perfectly-plastic  adhesive  shear  strain  distribution  in  a  stepped 
double-lap  joint  is  similar  to  solving  that  for  the  elastic  case  except  that  the  adhesive  under 
each  step  can  now  also  deform  fully  plastically  or  partly  plastically  and  partly  elastically. 


upper  bound 


lower  bound 


,  .  — ,1— -r. 


The  solution  assuming  perfectly-plastic  adhesive  deformation  on  the  kth  step  is 
as  given  by  equation  35; 
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7k(^]  =  -"20^  +C^+D.  (49) 

where  ^  is  a  lcx:al  coordinate  for  each  step.  The  boundary  condition  given  by 
equations  11  still  applies.  Hence,  on  the  first  step  the  constant  C  in  the  above 
equation  is: 


C=- 


P 

T|Ei 


(50) 


If  on  any  step  there  is  a  transition  from  plastic-to-elastic  adhesive  deformation 
then  the  shear  strain  and  its  derivative  given  by  equation  50  is  matched  with  the 
elastic  solution  (equation  37).  At  the  end  of  each  step  continuity  of  the  shear 
strain  is  enforced. 

The  derivative  of  the  shear  strain  at  the  start  of  the  kth  step  is  given  by 
equation  41 


nUEoto 


(51) 


The  traction  T„  at  the  start  of  kth  step  is  found  by  integrating  equation  1  for  the 
previous  k*l  steps  The  traction  T,  is  simply  equal  to  P-2T,.  Thus  once  the 
shear  strain  at  the  start  of  the  first  step  is  known  the  full  solution  can  be  found 

Once  again  a  computerised  iterative  scheme  was  implemented  whereby  repeated  guesses 
were  made  at  the  starting  shear  strain  until  the  load  transferred  was  equal  to 
p  E„t,. 

2  Eont„  +  E,  t, 

4.  COMPARISON  OF  THE  END-TAPERED  AND  THE  STEPPED  DOI  BLE-LAP 
JOINT  MODELS 

4.1.  Assuming  Elastic  Adhesive  Deformation. 

Table  6  contains  the  joint  parameters  used  in  the  model  of  an  end-tapered  double-lap  joint 
and  in  the  model  of  an  equivalent  stepped  double-lap  joint  assuming  elastic  adhesive 
deformation  The  adhesive  shear  strains  and  outer  adherend  strains  for  these  two  joints 
appear  in  the  figures  that  follow 
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Figures  15  and  16 

End-tapered  Double-Lap  Joint 

9  (degrees) 

2.7 

1  (mm) 

27 

(mm) 

0.13 

Stepped  Double-Lap  Joint 

to  (mm) 

0.13 

n  (no.  steps) 

10 

step  length  (mm) 

3 

Table  6. 


adhesive 

shear 

strain 


Figure  15. 
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Figure  16. 


Significantly,  however,  the  stepped  double-lap  joint  does  predict  a  lower  peak  shear  strain. 
The  outer  adherend  strains  are  close  in  an  average  sense. 

4.2.  Assuming  Elastic/Perfectly-Plastic  Adhesive  Deformation. 

Table  7  contains  the  joint  parameters  used  in  the  model  of  an  end-tapered  double-lap  joint 
and  in  the  model  of  an  equivalent  stepped  double-lap  joint  assuming  elastic/perfectly-plastic 
adhesive  deformation  The  adhesive  shear  strains  and  outer  adherend  strains  for  these  two 
joints  appear  in  the  figures  that  follow. 


Figures  17  and  18 

End-tapered  Double-Lap 
Joint 

0  (degrees) 

2.9 

1  (mm) 

27 

(™^i) 

0.15 

Stepped  Double-Lap  Joint 

to  (mm) 

0.15 

n  (no.  steps) 

10 

step  length  (mm) 

3 

Table  7. 


r 


24 

Allowing  the  adhesive  to  deform  in  an  elastic/perfectly-plastic  manner  leads  to  a  slight 
divergence  but  still  similar  prediction  of  the  adhesive  shear  strains  by  the  two  models.  The 
predicted  outer  adherend  strains  for  the  two  joints  remain  close  in  an  average  sense. 


adhesive 

shear 

strain 


Figure  1 7. 


$ 


Figure  16. 
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4.3.  Comparison  of  the  Two  Models  with  Measured  Outer-Adherend  Strains 

The  ability  of  the  two  models  to  predict  strain  in  the  outer  adherend  of  an  actual  joint 
(Figure  Ic)  was  investigated.  A  strip  of  strain  gauges  was  bonded  to  the  outer  adherend  of 
a  boron/epoxy  doubler  near  the  end  of  the  taper.  The  joint  was  loaded  and  the  strain 
measured.  The  results  and  the  predicted  strains  (calculated  using  the  elastic  models)  are 
shown  in  Figure  19. 


distance  along  joint  (mm) 


Figure  19. 

Clearly,  the  stepped  structure  of  the  doubler  is  reflected  in  the  strain  measured  on  the  outer 
ply.  The  stepped  double-lap  joint  model  is  better  suited  to  describing  the  lack  of 
smoothness  in  the  strain/distance  curve  though  the  end-tapered  model  is  accurate  in  an 
average  sense.  The  boron/epoxy  doubler  consisted  of  10  plies  and  a  step  length  of  3mm. 
The  resultant  taper  angle  was  about  2.5°.  Five  strain  gauges  were  bonded  to  the  doubler; 
the  strain  gauges  were  spaced  2mm  apart  starting  at  3mm  from  the  end  of  the  doubler. 

5.  CONCLUSION 

Two  alternative  mathematical  models  have  been  successfully  developed  to  predict  adhesive 
shear  strains  in  fibre-composite  reinforcements.  The  validity  of  both  the  end-tapered 
double-lap  joint  model  and  the  stepped  doubie-lap  joint  model  has  been  verified  by  an 
elastic  finite-element  analysis.  For  a  typical  repair  or  doubler  geometry  the  two  m^els 
yield  approximately  the  same  adhesive  shear  strains  for  both  elastic  and  elastic/perfectly 
plastic  analyses.  The  end-tapered  double-lap  joint  model  is  numerically  more  stable  than 
the  stepped  double-tap  joint  model  and  is  periiaps  therefore  more  useful  in  predicting 
adhesive  shear  strains.  The  stepped  double-lap  joint  model,  however,  gives  a  better 
prediction  of  the  strains  in  the  outer  adherend.  Estimation  of  the  peak  shear  strain  in  the 
stepped  double-lap  joint  may  be  simplified  by  umg  a  angle  ply  model  for  certain  joint 
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configurations.  The  stepped  double-lap  joint  geometry  appears  to  be  more  effective  for 
repair  than  the  end-tapered  double-lap  joint  geometry  as  the  predicted  peak  shear  stress  is 
lower.  The  two  models  could  be  useful  design  tools  when  used  in  conjunction  with 
adhesive  design  allowables  such  as  the  fatigue  threshold,  the  yield  strain,  and  the  strain  to 
failure. 
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It  ABSTRACT 

Bonded  fibre-composiie  doublers  are  increasingly  being  used  to  reinforce  and/or  repair  damaged  or 
underdesigned  metallic  aircrerft  structure.  This  trend  will  continue  as  the  average  age  of  Australian  dvilian 
and  military  aircraft  increases.  Two  mathematical  models  of  the  stress  state  in  bonded  doublers  are 
presented  in  this  report:  the  end-tapered  double-lap  joint  and  the  stepped  double-lap  joint.  The  mathematical 
develoi  ’ment  of  the  two  models  is  detailed  and  the  fidelity  of  the  predicted  stress  states  compared  with  that 
obtained  from  finite  element  analyses  and  experiment. 
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